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Equilibrium properties of the three-dimensional isotropic Heisenberg spin glass are studied by 
extensive Monte Carlo simulations, with particular attention to the nature of its phase transition. 
A finite-size-scaling analysis is performed both for the spin-glass (SG) and the chiral-glass (CG) 
orders. Our results suggest that the model exhibits the CG long-range order at finite temperatures 
without accompanying the conventional SG long-range order, in contrast to some of the recent works 
claiming the simultaneous SG and CG transition. Typical length and time scales which represent 
a crossover from the spin-chirality coupling regime at short scales to the spin-chirality decoupling 
regime at long scales are introduced and examined in order to observe the true asymptotic transition 
behavior. On the basis of these crossover scales, discussion is given concerning the cause of the 
discrepancy between our present result and those of other recent numerical works. 
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I. INTRODUCTION 

Spin glasses (SGs) have attracted the attention of re- 
searchers both in experiments and theory as a prototype 
of complex systems with quenched randomness i SGs are 
random magnets in which magnetic ions interact with 
each other either ferromagnctically or antiferromagneti- 
cally, depending on their positions. Most of theoretical 
works have been so far devoted to the minimal SG model, 
i.e., the Ising Edwards and Anderson (EA) model. After 
a discussion early in the 1980s, it is now widely believed 
that a three-dimensional (3D) Ising SG model exhibits a 
SG phase transition at a finite temperature. Large-scale 
Monte Carlo (MC) simulations presented an evidence for 
the finite-temperature SG ordcring^^ Subsequently, the 
critical exponents evaluated by MC simulations were con- 
sistently compared with those evaluated experimentally 
for the Ising-like SG compound FeMnTiOs. 

Compared to the Ising case, the nature of the phase 
transition of continuous spin systems such as XY and 
Heisenberg SGs arc still poorly understood. Since many 
SG magnets including canonical SG possess only weak 
magnetic anisotropy, an isotropic Heisenberg SG model, 
rather than the strongly anisotropic Ising model, is ex- 
pected to be a realistic model of SG magnets. Exper- 
imentally, an equilibrium SG phase transition has been 
established in real SG materials via measurements of the 
divergent non-linear susceptibility, etc. In sharp con- 
trast to experiments, earlier theoretical studies on the 
Heisenberg SG model indicated that the standard SG 
long-range order occurred only at zero temperature in 
three dimensions4*5i£iL£ 

In order to solve this apparent puzzle, a chirality mech- 
anism of experimentally observed SG transitions was pro- 
posed by Kawamura&i^. This scenario is based on the 
assumption that an isotropic 3D Heisenberg SG exhibits 
a finite-temperature chiral-glass (CG) transition without 



the conventional SG order. In terms of symmetry, among 
the global symmetries of the isotropic Hamiltonian, only 
the Z2 spin-reflection (or spin-inversion) symmetry as- 
sociated with the chirality is spontaneously broken with 
keeping the 50(3) spin-rotation symmetry Indeed, some 
numerical studiesAi*i£ claimed that the standard SG or- 
der associated with the freezing of the Heisenberg spin 
occurred at a temperature lower than the CG transition 
temperature, i.e., Tsq < Tqg, possibly with TgG = 0. 
It means that the spin and the chirality are decoupled 
on long length and time scales, although the chirality is 
locally defined as a composite operator of the spin vari- 
ables. 

In this chirality scenario of experimental SG transi- 
tions, essential features of many of the real SG tran- 
sition and of the SG ordered state are determined by 
the properties of the CG transition and of the CG state 
of the fully isotropic system. The role of the magnetic 
anisotropy is secondary which re- couples the spin to the 
chirality and reveals the CG transition in the chiral sector 
as an anomaly in experimentally accessible spin-related 
quantities. The scenario successfully explained the phase 
diagram under magnetic fields observed by the recent nu- 
merical simulationiii^ and experiments^ 

More recently, however, some researchers argued a pos- 
sibility that in the 3D Heisenberg SG model the spin 
ordered at a finite temperature simultaneously with the 
chirality, i.e., T SG = T CG > Q, 16 i 17 i 18 i 19 i 20 Thus, the na- 
ture of the ordering of the 3D Heisenberg SG, as well as 
the validity of the chirality scenario, is now under de- 
bate. Under such circumstances, it is highly interesting 
to perform further extensive numerical studies of the 3D 
Heisenberg SG in order to clarify the true nature of its or- 
dering. In the present study, we investigate both the SG 
and CG orderings of the model by means of a large-scale 
equilibrium MC simulation. 

Interestingly, recent experiments reported on a quali- 
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tative difference in aging phenomena between a canonical 
Heisenberg-like SG and an Ising-like SGi^i We also ex- 
pect that the full understanding of the equilibrium prop- 
erties of the 3D Heisenberg SG will also give a valuable 
insight into these off-equilibrium properties of SGs. 

The article is organized as follows. In Sec.[nj we give a 
background of the present numerical study. In Sec. Ill Al 
we explain first the basics of the chirality mechanism. 
In Sec. IIIBI wc introduce the crossover length and time 
scales beyond which the spin and the chirality are decou- 
pled with each other. These length and time scales are 
crucially important in the chirality mechanism, and are 
also essential in properly interpreting the numerical data 
of MC simulations. In Sec. IIIII we explain the model 
and the MC method employed. In Sec. II VI we introduce 
various physical quantities measured in our MC simula- 
tions, while the results of our simulations are presented 
in Sec. In view of our MC results, we examine and 
discuss in Sec. lVII the recent numerical results on the 3D 
Heisenberg SG by other authors. Finally, we present a 
brief summary of the results in section VII. 



II. BACKGROUND 

In this section, we wish to give a background of the 
present numerical study of the 3D isotropic Heisenberg 
SG. First, we explain the basics of the chirality mecha- 
nism of experimental SG transition as proposed in Refs.0 
and 03 Then, we explain the notion of the spin-chirality 
decoupling, together with the crossover length and time 
scales which play a crucially important role in the chi- 
rality mechanism and are also essential in properly inter- 
preting the numerical data. 



A. Chirality mechanism 

Chirality is an Ising-like multi-spin variable represent- 
ing the sense or the handedness of noncollinear spin struc- 
tures induced by spin frustration. In frustrated magnets 
with continuous spins, the chirality often plays an essen- 
tial role in their magnetic ordering. The local chirality 
Xip, at the ith site in the /x-direction may be defined by 

to^S^/fSjxS,^), (1) 

e^O-* = x i Vi z ) being a unit lattice vector along the fi axis. 
This quantity is often called a scalar chirality: It takes 
a non-zero value only when the three neighboring spins 
take the non-coplanar configuration in spin space, while 
it vanishes for the collinear or the coplanar spin config- 
uration. The chirality defined above is a pseudoscalar 
variable since it is invariant under the global SO (3) spin 
rotations but changes its sign under the global Zi spin 
reflections or inversions. 



The chirality mechanism of Refs.l^andlTol takes the fol- 
lowing two-step strategy in explaining the real SG transi- 
tion: The first step concerns with the property of the fully 
isotropic Heisenberg SG, an idealization of experimental 
SG materials. The chirality scenario claims that the fully 
isotropic Heisenberg SG exhibits a finite-temperature CG 
transition without the conventional SG long-range or- 
der. The CG transition breaks only the Z2 spin-reflection 
symmetry with keeping the 5*0(3) spin-rotational sym- 
metry. The occurrence of the CG transition necessarily 
entails the spin-chirality decoupling. 

Obviously, such a scenario does not apply to the infi- 
nite dimensional limit, i.e., to the mean-field Heisenberg 
Sherrington-Kirkpatrick (SK) model, in which the spin 
itself, not the chirality, behaves as an order parameter of 
the transition. Due to the noncoplanar nature of the spin 
configuration in the SG state, the SG long-range order 
trivially accompanies the CG long-rage order, whereas 
the opposite is not necessarily true. One should note 
that, in the conventional case where the spin variable 
is a proper order parameter of the transition as in the 
case of the SK model, the chirality, which is given by 
the multiple of the spin, exhibits a less singular behavior 
than the spin at the SG transition. In fact, the chirality 
shows only moderate behavior at the SG transition of the 
mean-field Heisenberg SK model in which the spin, not 
the chirality, is the order parameter of the transition,^ 

In contrast to the mean-field model or the high- 
dimensional Heisenberg SG models, the problem could 
be very non-trivial in lower dimensions where the order- 
parameter fluctuation might change the nature of order- 
ing dramatically. At present, there seems to be no con- 
sensus about the lower critical dimension dl£P of the SG 
order, while the corresponding upper critical dimension is 
expected to lie six. The CG order, if exists, may emerge 
slightly above, at, or below dg£ D . 

It has been proved that the SG long-range order 
does not exist at any finite temperature in the two di- 
mensional Heisenberg SG^ The numerical domain-wall 
renormalization-group calculation as well as the MC sim- 
ulation suggested that both the spin and the chirality or- 
dered only at zero temperature^. Interestingly, however, 
the estimated SG and CG correlation-length exponents 
at this T = transition differ significantly from each 
other, i.e., i/qg > ^sg^* This implies that in 2D the 
spin and the chirality are decoupled at long length scale, 
the chirality dominating the long-length behavior. 

In view of these transition behaviors of the Heisenberg 
SGs, it appears likely that the principal player in long- 
scale phenomena changes from the spin to the chirality 
as the spacial dimensionality is decreased. Thus, the be- 
havior in dimension three is the current issue, which is 
the subject of the present work. 

The second step of the chirality mechanism concerns 
with the effect of the random magnetic anisotropy which 
inevitably exists in real SG magnets. The random 
anisotropy energetically breaks the SO(3) spin-rotation 
symmetry in the Hamiltonian, retaining the Z2 inversion 
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FIG. 1: A schematic figure of the crossover between the spin- 
glass and chiral-glass correlation lengths (correlation times) 
expected in the chirality mechanism. According to the chiral- 
ity mechanism, the CG correlation length (correlation time) 
diverges toward the CG transition temperature T — Tcg, 
while the SG one diverges at a lower temperature T = Tsg < 
Tcg- 



symmetry only. When the anisotropic system exhibits 
the CG long-range order with spontaneously breaking 
the Z2 inversion symmetry, there no longer remains any 
global symmetry degree of freedom to leave the system in 
the paramagnetic phase. Hence, once the CG order oc- 
curs in the presence of the random anisotropy, the spin 
degree of freedom also behaves like the chirality. This is 
the spin-chirality rccoupling due to the random magnetic 
anisotropy. 

Such an anomaly revealed in the spin sector via the 
random magnetic anisotropy can be detected experimen- 
tally by standard magnetic measurements e.g. , as a di- 
vergence of the nonlinear susceptibility etc , whereas the 
CG long-range order is difficult to observe experimen- 
tally. 

We note that in this mechanism the anisotropy plays 
only a secondary role: The anisotropy certainly reduces 
the symmetry of the Hamiltonian relative to the fully 
isotropic system, but does not change the broken sym- 
metry of the transition. The critical properties of the 
CG transition and of the low-temperature CG phase are 
expected to be not affected by the magnetic anisotropy, 
which, however, are now directly observable via the stan- 
dard spin-related quantities. This chirality scenario pre- 
dicts that experimentally observed SG transitions belong 
to the same universality class as that of the the CG tran- 
sition of the fully isotropic model. It is thus highly inter- 
esting to clarify the nature of the phase transition of the 
ideal isotropic Hciscnbcrg SG. 



B. Spin-chirality decoupling/coupling scenario 

As mentioned above, in lower dimensions, a relevant 
degree of freedom which dominates the long-scale phe- 
nomena might well change from the spin to the chirality. 
The chirality scenario expects that in 3D there exists a 
crossover temperature T x which separates the two tem- 
perature regimes, as illustrated in Fig. ^ I n the higher 
temperature regime, the SG correlation length is longer 
than the CG correlation length, dominating the long- 
scale phenomena. This is simply due to the fact that 
the sensible definition of the local chirality requires the 
development of the spin short-range order of at least a 
few lattice spacings. As the temperature is decreased, 
both the SG and CG correlation lengths grow, but at 
different rates, so that the CG correlation length eventu- 
ally outgrows the SG correlation length at the crossover 
temperature T x . An example of such a crossover behav- 
ior between the spin and the chiral correlations can be 
seen exp licitly in a certain toy model: See Fig. 10 of 
Ref. Then, the relevant degree of freedom for the 
long-scale behavior changes at T x from the spin to the 
chirality. Below T x , the long-scale phenomena are gov- 
erned by the CG correlation, not by the SG one. This 
is the spin-chirality decoupling expected to occur in the 
fully isotropic model. 

Let us discuss in some detail the finite-size effect in- 
herent to the simulation data in the critical region. The 
situation here is not simple because the system has two 
length scales, each associated with the spin and with 
the chirality. Suppose that the CG transition occurs at 
T = Tqg without the conventional SG long-range or- 
der. Then, the crossover temperature T x at which the 
CG correlation length outgrows the SG correlation length 
should be located somewhat above Tcg: See Fig.l. A 
necessary condition for detecting the spin-chirality de- 
coupling is that the measurement temperature lies be- 
low T x . It is, however, not enough. At a temperature 
below T x , one needs to probe the system beyond the 
crossover length above which the spin-chirality decou- 
pling becomes eminent. Thus, a large-size simulation 
exceeding the crossover length is required in order to 
detect the spin-chirality decoupling. Unfortunately, the 
crossover length scale is unknown a priori, and is to be in- 
vestigated by numerical simulations. A natural criterion 
might be that it is given by the SG correlation length at 
the crossover temperature T x as shown in Fig. ^ Even 
in the CG ordered phase, the spin-chirality decoupling 
might hardly be observable at the length scale below the 
crossover length. Rather, it is natural to expect that 
the trivial spin-chirality coupling is observed below the 
crossover length scale because the chirality is a composite 
operator of the spin on the short scale of lattice spacing 
not independent of the spin, roughly being x ~ S 3 . 

In the CG critical region, the chirality-related quanti- 
ties should exhibit the true asymptotic critical behavior, 
e.g., a power-law singularity characterized by the asso- 
ciated CG exponents. At short length scales below the 
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FIG. 2: A schematic figure of the time evolution of the spin- 
glass and chiral-glass two-time autocorrelation functions. 



crossover length scale, due to the trivial coupling between 
the spin and the chirality, even the spin-related quanti- 
ties are expected to exhibit the similar critical behavior 
to the chirality-related quantities. Namely, up to the 
crossover length scale, it seems as if the SG order devel- 
oped as a long-range order. It is intrinsically difficult at 
shorter length scales to distinguish such a pseudo-critical 
behavior induced by the CG long-range order from the 
true SG long-range order. Hence, it is crucially impor- 
tant to estimate the crossover length scale and to study 
the long-scale behavior of the system beyond this length 
scale. 

Essentially the same argument also applies to the case 
of the temporal scale. As an example, we discuss here 
the behavior of the autocorrelation functions based on 
the notion of the crossover time scale. Figure shows 
the schematic representation of the behavior of the SG 
and CG autocorrelation functions below Tqg expected 
from the spin-chirality coupling/decoupling picture. In 
the long-time limit, the chirality autocorrelation func- 
tion is expected to saturate to a certain finite value after 
an initial fast decay, while the spin autocorrelation func- 
tion is expected to decay toward zero asymptotically. A 
comment is in order concerning the transient behavior of 
the spin autocorrelation function at time scales shorter 
than the crossover time scale: At shorter time scales, the 
spin autocorrelation function might exhibit the pseudo- 
ordering feature dictated by the CG one through the 
trivial spin-chirality coupling. This might well lead to a 
hump-like pseudo-ordering structure in the time depen- 
dence of the spin autocorrelation function as shown in 
Fig. 2, which, however, does not persist in the long-time 
limit beyond the crossover time scale. This means that, 
in order to properly discuss the true asymptotic behav- 
ior of the dynamics of the model, a particular attention 
should be paid to the crossover time scale. 



III. THE MODEL AND THE MONTE CARLO 
METHOD 

We study a classical Heisenberg model defined by the 
Hamiltonian, 

H(S) = -^2JiA-S j , (2) 

m 

where Si = (Sf, , Sf) is a three-component unit vector, 
and the summation runs over all nearest-neighbor pairs. 
The lattice is a simple-cubic lattice with the total number 
of N = L 3 sites. The nearest-neighbor couplings take 
the values ± J randomly with equal probability. Periodic 
boundary conditions are imposed for all the directions. 
The lattice sizes studied are L = 8, 12, 16 and 20, where 
the sample average is taken over 976 (L = 8), 964 (L 
12:. 280 (L = 16), and 32 (L = 20) independent bond 
realizations. 

We perform an equilibrium MC simulation of the 
model. In our simulation, we make use of the exchange 
MC method^ which is also called parallel tempering^ 
In the exchange MC method, one MC step consists of 
two elementary updates, a standard single-spin heat-bath 
flip£ and an exchange trial of spin configurations at neigh- 
boring temperatures. The latter reduces the slow re- 
laxation at low temperatures with the help of the high- 
temperature fast dynamics. The method has turned out 
to be quit efficient in thcrmalizing a wide class of hardly- 
relaxing systems such as SG systems and proteins. We 
ensure equilibration by checking that various observables 
attain stable values, no long er changing with the amount 
of MC steps: See Refs. if and |H for further details of 
the equilibration procedure. Our MC simulations have 
been performed up to the size L = 20 and up to the tem- 
perature T/J — 0.15. This could be achieved only by 
using the exchange MC method. The numbers of tem- 
perature points used in our exchange MC method are 32 
for L = 8, 12, 16, and 48 for L = 20. 

Error bars are estimated via sample-to-samplc fluctu- 
ations for the linear quantities such as the order parame- 
ters, and by the jackknife method for the non-linear quan- 
tities such as the Binder parameter and the correlation 
length mentioned below. 



IV. PHYSICAL OBSERVABLES 

In the present section, we introduce various physical 
quantities observed in our simulations, and discuss some 
of their basic properties. 

In glassy systems, it is often convenient to define as 
an order parameter an overlap variable between two in- 
dependent systems with the same Hamiltonian. For the 
Heisenberg spin, the overlap may be defined as a tensor 
variable between the y, and v components (p, v = x,y, z) 
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of the Heisenberg spin by 



while that for the CG oder is defined by 



1 N 



(3) 



ffCG 



(ii) 



where the upper suffixes (1) and (2) denote the two repli- 
cas of the system with the same interaction set. 

A chiral overlap is defined in terms of the local chiral 
variable (|TJ by 



Qx 



— Vy (1) y (2) 



(4) 



The squared SG order parameter is then given by 



(2) 
%G = 
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where (■ ■ • ) denotes a thermal average and [• • • ] denotes 
an average over the bond disorder. The corresponding 
squared CG order parameter is defined by 



(2) 

9cg 
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which is normalized by the mean-square amplitude of the 
local chirality, 
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In the thermodynamic limit, these Binder parameters 
are normalized to unity in the non-degenerate ordered 
state, and to zero in the high-temperature disordered 
state. Since the Binder parameter is a dimensionless 
quantity, and the dimensionless quantity should be size- 
independent at the critical temperature T c , the Binder 
parameters of different system sizes plotted as a function 
of temperature should yield a crossing or merging point 
at T c . 

In terms of the fc-dependent overlap variable, one can 
define the Fourier-transformed two-point CG and SG cor- 
relation functions. For the CG, the fc-dependent chiral- 
overlap is defined by 



N 



Qx( k ) 



^Ex^X^exp^fc-n), (12) 



N 
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in which the chiral variable along the cc-axis, Xixi is con- 
sidered. The Fourier-transformed CG correlation func- 
tion is then defined by 



&)= (l<fc(*')l 



(13) 



For the SG, the fc-dependent spin-overlap is defined by 



The local chirality amplitude remains non-zero only when 
the spins have a non-coplanar structure locally. This 
quantity weakly depends on the temperature, in contrast 
to the Heisenberg spin variable whose amplitude is fixed 
to be unity by definition. In the high-temperature sym- 
metric phase, these SG and CG order parameters are 
essentially equivalent to the associated SG and CG sus 
ccptibilities defined by xsg 
respectively. 

A standard finite-size scaling of the second-order tran- 
sition for the equilibrium SG and CG order parameters 
takes the form 



NQsg and Xcg 



.,(2) 
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-(i+v) 



'( 



T-TJL 



l/u 
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where v is the exponent of the correlation length, and r\ 
is the exponent describing the decay of the correlation 
function at the critical point T — T c . At T c , the order 
parameter decays as a power law with the size L, 



q (2) K £-(!+•»>. 



(9) 



One often uses the Binder parameter to estimate the 
critical temperature. In the Heisenberg SG, the Binder 
parameter for the SG order is defined by 



i/ t1 Q K<z 4 >] 
9SG = 2{ U - 9 WW 



q 2 = 
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whereas the Fourier-transformed SG correlation function 
is defined by 



E 



(15) 



Via these CG and SG Fourier-transformed correlation 
functions, the associated CG and SG finite-system corre- 
lation lengths are defined by 



i / gW(o) 

2sin(fc m /2)Y 5 (2)(fc m ) 



(16) 



where fc m = (2tt/L,0, 0) and fc m = |fc m |. 

One can then define a dimensionless quantity, the nor- 
malized correlation length Z,cg/L and £sg/L. Since £/L 
is dimensionless, it should exhibit the same scaling prop- 
erty as the Binder parameter. Thus, the ratio £/L for 
different L should cross or merge at the critical temper- 
ature. 

In probing the nature of the low-temperature glassy 
ordered phase, one useful quantity is the distribution of 
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the overlap. The chiral-ovcrlap distribution is defined by 



(17) 



(2) 

The squared CG order parameter defined above 
is the second moment of the chiral-overlap distribution 
function. 

The spin-overlap distribution is defined originally in 
the tensor space with 3x3 = 9 components. To make 
the quantity more easily tractable, one may define the 
diagonal spin overlap which is the trace of the original 
tensor overlap, and introduce the associated diagonal- 
spin-ovcrlap distribution by 



- E 




(18) 



This distribution function is symmetric with respect to 
<Zdiag = 0, and is expected to be a Gaussian distribu- 
tion around (faiag = in the high-temperature disor- 
dered phase. Reflecting the fact that the diagonal-spin- 
overlap is not invariant under the global 0(3) spin rota- 
tion, P(<7diag) in the possible SG ordered phase develops a 
nontrivial shape, not just consisting of the delta-function 
peaks related to q^A, even when the ordered state is 
a trivial one simply described by a self-overlap qEA~ 
If the possible SG ordered state accompanies a replica- 
symmetry breaking (RSB), further nontrivial structures 
would be added to -P(<Zdiag)- Meanwhile, it is recently 
shown in Ref. 1 221 that, in the possible SG ordered state, 
the diverging peaks corresponding to the self-overlap nec- 
essarily appears in -P(gdiag) at gaiag = ±g<ZEA in the ther- 
modynamic limit. Hence, the existence or nonexistence 
of these divergent peaks could be used as an unambigu- 
ous measure of the possible SG long-range order in the 
Heisenbcrg SG, irrespective the occurrence of the RSB. 

Another interesting feature of glassy systems might be 
the nature of their sample-to-sample fluctuations, partic- 
ularly its possible non-self-averageness. As an indicator 
of the lack of self-averageness, one may use the so-called 
A parameter^ For the CG order, it is defined by 



A CG (T) 
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while for the SG order, 
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The order parameter is non-self-averaging when the as- 
sociated A parameter is non-zero, and is self-averaging 
when A is equal to zero. 

One can also define the so-called Guerra parameter 
G~ For the CG order, it is defined by 



Gcg(T) 



Kg 2 ) 2 ] 



(21) 



while for the SG order, 
Gsg(T) : 



[(g 2 ) 2 ] - [(g 2 )] 2 

K<7 4 >] - Kg 2 )] 5 



(22) 



Unlike the A parameter, the G parameter can take a 
non-zero value even when the ordered state is a trivial 
one without accompanying the RSB. 31-32 The G param- 
eters are related to the A parameters and the Binder 
parameters g via the relations, 
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These relations indicate that, so long as the Binder pa- 
rameter g takes any value different from unity in the or- 
dered phase, a non-zero A necessarily means a non-zero 
G. By contrast, if the Binder parameter g takes a value 
unity in the ordered phase, a non-zero A may or may not 
mean a non-zero G. 

Information about the equilibrium dynamics can be 
obtained from the spin and chiral two-time autocorrela- 
tion functions defined by 



C s {t) 



c x {t) = 



N ^ 

-E 



+ to))], 



(25) 



[(XiMxi»(t + to))}, (26) 



where the time evolution in our MC simulation is made 
according to the standard heat-bath updating not accom- 
panying the exchange process. Initial spin configurations 
at t = t$ are taken from equilibrium spin configurations 
generated in our exchange MC runs. Below the transition 
temperature T c (if any), these autocorrelation functions 
converge in the long time limit to the Edwards- Anderson 
SG and CG order parameters, whereas above T c these au- 
tocorrelation functions decay exponentially toward zero 
with a characteristic correlation time, which diverges as 
the temperature T approaches T c . Just at T c , the auto- 
correlation functions exhibit a power-law decay, 

t -/3/zu 



C{t) 



(27) 



where z is the dynamical critical exponent. These fea 
tures are described by the standard bulk dynamical seal 
ing form, 

\T-T c ff(t\T-T c D, 



C(t) 



(28) 



where f(x) is a scaling function whose asymptotic forms 
for x 1 and x 3> 1 are x~^l zv and exp(— a;), respec- 
tively. 



V. THE NUMERICAL RESULTS 

In the present section, we show our MC results of the 
three-dimensional isotropic Heisenbcrg SG model. 
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FIG. 3: The temperature and size dependence of the chiral- 
glass Binder parameter; upper figure (a), and of the spin-glass 
Binder parameter; lower figure (b), of the 3D ±J Heisenberg 
SG. 



A. Binder parameter 

In Fig. 13 we show the temperature and size depen- 
dence of the Binder parameters both for the chirality; 
upper figure (a), and for the spin; lower figure (b). As 
can be seen from Fig.[3{a), a crossing of the CG Binder 
parameter gcG of different L is observed on the negative 
side of gcG, not on the positive side as in the standard 
cases. With increasing L, the crossing temperature grad- 
ually shifts toward lower temperatures. A behavior sim- 
ilar to this has also been observed in the Binder param- 
eter of other models, including the Heisenberg SGiS^ 
and the mean-field SG2&24. In particular, in a class of 
mean- field SG models exhibiting a one-step RSB, the 
Binder parameter at the transition point T c takes a neg- 
ative value, sometimes even negatively divergent ^2*24 It 
implies that the temperature at which the Binder pa- 
rameter for finite L takes a minimum, a dip tempera- 
ture Tdip(.L), approaches the critical temperature, i.e., 
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Tdip{L) — > T c as L — > oo. Recently, this method of 
estimating the bulk transition temperature was success- 
fully applied to the Heisenberg SGii^S In Fig.0] we plot 



FIG. 4: The dip temperature of the chiral-glass Binder pa- 
rameter gcG is plotted against l/L. The solid line represents 
a linear fit of the data. Its extrapolation to the L — > oo limit 
gives an estimate of the bulk chiral-glass transition tempera- 
ture, Tcg/J = 0.194(5). 



Tdip(L) against l/L. An extrapolation to the thermody- 
namic limit l/L — > gives us an estimate of the bulk CG 
transition temperature, Tcg/ J = 0.194(5). 

By contrast, as shown in Fig.^Jb), the SG Binder pa- 
rameter gsG monotonically decreases toward zero with 
increasing L at all temperatures studied. There is no 
signature of the transition in the investigated tempera- 
ture range, no negative dip nor the crossing, in contrast 
to the CG Binder parameter. This suggests that the SG 
transition temperature, if any, is located at a temperature 
lower than the temperature range studied here. Fig.^b) 
reveals, however, that an anomalous bend appears in gsc 
for larger sizes L > 16 at around T/J ~ 0.22, close to the 
CG transition temperature, although <?sg never becomes 
size-invariant at any temperature, as it should have been 
in a second-order transition. The reason why gsG ex- 
hibits such an anomalous bend around Tcg might be 
understood as follows: At the CG transition, a reflection 
symmetry is spontaneously broken and the entire phase 
space is divided into ergodic components, in each of which 
a proper-rotational symmetry is still preserved. As a re- 
sult, the ordering behavior of the Heisenberg spin would 
change at Tcg > though the spin itself does not order even 
below Tcg- We note that a similar bend in gsG has also 
been observed in the two-dimensional Heisenberg SG™ 
where the absence of a finite-temperature SG transition 
has been well established™ 

It is sometimes argued in the literature that the 
Binder-parameter analysis might not work in the SG 
problem. Such a suspicion might partly be based on 
the observation that only weak merging behavior was 
observed at the SG transition temperature of the three- 
dimensional EA Ising model which is believed to exhibit 
a finite-temperature SG transitioni2ii2^ As long as the 
SG long-rage order really sets in at finite temperatures, 
however, it is hardly conceivable that the Binder pa- 
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rameter for asymptotically large lattices exhibits a non- 
singular behavior only. In particular, the Binder param- 
eter should become scale-invariant at the SG transition 
point, so long as the transition is continuous. Indeed, in a 
recent MC simulation of the mean-field Heisenberg SG^2. 
which is known to exhibit a non-zero SG long-range order 
below TsGi a clear crossing of the SG Binder parameter 
gsG has been observed at Tsg, in sharp contrast to our 
present data of Fig.|3fb). 

Then, one might argue that the finite-size effect would 
be significant here in ggQ and the largc-T asymptote 
might still be far away. One sees from Fig.^a), however, 
that the CG Binder parameter gcG for our two largest 
sizes L = 16 and 20 exhibits an almost size-invariant 
behavior at and below Tcg- If the Heisenberg spin or- 
ders simultaneously with the chirality, and if the spin is 
the order parameter of the transition and the chirality is 
only composite (secondary), it seems a bit hard to under- 
stand why the chirality exhibits an almost scale- invariant 
near-critical behavior for L > 16, while the Heisenberg 
spin still exhibits an off-critical scale-dependent behavior. 
Hence, the behavior of g^G observed in Fig.^b) remains 
to be resolved if the occurrence of a finite-temperature 
SG transition is to be accepted in the investigated tem- 
perature range. 



B. Order parameter 

In Fig. [3 we show the size dependence of the squared 
CG and SG order parameters, q^Q and q^Q, for vari- 
ous temperatures. In the upper figure (a), a double- 

(2) 

logarithmic plot of the CG order parameter <j CG is shown 
against the system size L. One generally expects that at 

(2) 

T c the data of q C Q should lie on a straight line. In fact, 

(2) 

the data of <7 CG show a clear straight-line behavior around 
T/J = 0.19, which is close to the CG transition temper- 
ature obtained by our analysis of gec- The critical-point 
decay exponent ?7cg can be estimated from the slope of 
this straight line, yielding 1 + ?7cg ~ 1-8. At higher tem- 
peratures, a deviation from the straight line, a downward 
trend, is observed indicative of the disordered phase. At 
lower temperatures, particularly at our lowest tempcra- 

(2) 

ture simulated T/J = 0.15, the data of g CG show a clear 
upward trend. This suggests that this temperature is in- 
deed below Tcg, and that the low-temperature phase is 

(2) 

a rigid one characterized by a non-zero <? C G> n °t likely 
to be a critical phase like the Kosterlitz-Thoulcss (KT) 
phase. 

In Fig. EJb), a double-logarithmic plot of the corre- 

(2) 

sponding SG order parameter q^ is shown against the 
system size L. Again, the data are expected to lie on a 
straight line at the critical SG transition temperature, if 
any. Such a straight-line behavior, however, is found only 
at our lowest temperature studied T/J = 0.15, whereas 

(2) 

%G never exhibits an upward trend characteristic of the 
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FIG. 5: Double-logarithmic plot of the squared chiral-glass 
order parameter Qqq; upper figure (a), and of the squared 
spin-glass order parameter gi^; lower figure (b), as a function 
of the system size L for several temperatures around the ex- 
pected chiral-glass transition temperature. Straight lines are 
drawn by connecting the two data points of L = 8 and 12 at 
each temperature. 



long-range ordered phase at any temperature studied, in 

(2) 

sharp contrast to the behavior of q CG - At temperatures 
higher than T/J = 0.15, including the one at the CG 

(2) 

transition temperature T/J ~ 0.19, the data of q§Q 
show a linear behavior for smaller sizes, which gradu- 
ally changes into a downward trend for larger sizes. This 
can simply be interpreted as a size-crossover which oc- 
curs around the length scale of the SG correlation length 
at each temperature. We note that such a size-crossover 
is clearly discernible even at a temperature T/J = 0.17 
which is below Tcg- The length scale of the crossover, 
comparable to the spin correlation length, grows as T 
decreases, and it is considered to exceed our largest size 
L = 20 at around T/J = 0.15. 

This observation strongly suggests that the standard 
SG transition temperature of the model is lower than 
T/J = 0.15, and that, at least within the temperature 
range 0.15 < T < 0.19, solely the CG long-range order 
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FIG. 6: The curvatures of the L-dependence of the squared 
spin-glass order parameter </gQ (open symbols), and of the 
squared chiral-glass order parameter q^Q (filled symbols) , are 
plotted versus the temperature. The curvature is expected to 
be zero at the respective transition temperature. 

exists without the standard SG long-range order, i.e., one 
has T CG > T SG . 

To make the situation more pronounced, we esti- 
mate following Refs. and the curvatures of the L- 

(2) (2) 

dependence of the two order parameters, q§Q and <7 CG , 
via second-order polynomial fits to the data of Figs. [S] 
The curvature is expected to be zero at the respective 
transition temperature. As shown in Fig. El the curva- 
ture for the CG crosses the zero-axis around T/J ~ 0.19, 
while that for the SG does not cross the zero-axis there, 
but marginally touches on it at a lower temperature, 
T/J ~ 0.15. The result indicates that the two transi- 
tion temperatures, Tsg and Tqq, are well separated. 



C. Finite-size scaling of the order parameter 

In order to estimate the correlation-length exponent 
v associated with the CG transition, we apply the stan- 
dard finite-size scaling analysis to the squared CG or- 
der parameter g CG based on Eq. ©. By taking \T — 
TcgI/TcgL 1 /" as the scaling variable, the best data col- 
lapse is obtained with Tqg/J = 0.19, ^cg = 1-2 and 
?7cg = 0.8. As shown in Fig. H the data both below and 
above T C G scale fairly well. If |1/T - l/T CG |?CG£ 1/iy 
is taken as the scaling variable, on the other hand, a 
slightly larger value of v, i.e., vqg = 1.4 and r]cG = 0.7, 
is preferred. The observed difference in the best values of 
the exponents might be due to the correction to scaling. 
Thus, we finally quote T CG /J = 0.19(1), v CG = 1.3(2) 
and t?cg = 0.8(2). The error bar is estimated by ex- 
amining the quality of the fits with varying the scal- 
ing parameters. The estimated values of critical expo- 
nents are compatible with the previous values obtained 
for the Hciscnberg SG model but with the Gaussian bond 



FIG. 7: Finite-size-scaling plot of the squared chiral-glass 
order parameter. The best scaling is obtained with Tcg/J = 
0.19, v = 1.2 and 1 + 77 = 1.8. 



distribution^ and also with those for the ± J Heisenberg 
SG under external field o 1 ?'?? . 

After establishing the occurrence of a finite- 
temperature CG transition, we next wish to re-examine 
via the finite-size scaling analysis the issue whether the 
standard SG order occurs at the same temperature with 

the CG order or not. In Fig. [HI we show a finite-size 

(2) 

scaling plot of the SG order parameter q§Q, assuming 
a simultaneous CG and SG transition with a common 
correlation length exponent, i.e., we set T§q/J = 0.19 
and z^sg = 1.2. Although the data turn out to scale well 
for smaller sizes, a significant deviation from the scaling 
is seen for larger sizes and at lower temperatures. The 
quality of the scaling is not improved if one tries to ad- 
just v to somewhat larger values. A similar poor scaling 
behavior is also observed even when one instead chooses 
|1/T — I/TqgITcgL 1 ^ as the scaling variable, and tries 
to adjust the scaling parameters around vsq = 1.4. The 
data for smaller sizes turn out to scale best with choos- 
ing rysc = —0.1. These parameter values ^sg = 1.2 and 
Vsg = — 0.1 are close to the values reported by Matsub- 
ara et al in Ref. l39l Hence, for the SG order parameter, 
we have observed a pseudo-critical behavior for smaller 
sizes, as well as a systematic deviation from the scaling 
for larger sizes. If the observed deviations were due to 
the correction-to-scaling, the scaling should be better for 
larger sizes, which is opposite to our present observa- 
tion. Therefore, we do not consider the apparent scaling 
obtained for smaller sizes with Tsg/ J = 0.19 to be ac- 
ceptable as a true asymptotic scaling. 

In Fig. |2| we show a finite-size-scaling plot of the SG 
order parameter q^l using the same data as in Fig. [H] 
but now assuming a zero-temperature SG transition, i.e., 
Tsg = and t]sg = — 1- The value r\ = — 1 is generically 
expected for a zero-temperature transition with the non- 
degenerate ground state. As shown in Fig. [51 the best 
data collapse is obtained by choosing v^g = 2.2. If one 
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FIG. 8: Finite-size-scaling plot of the squared spin-glass or- 
der parameter, assuming a simultaneous chiral-glass and spin- 
glass transition with a common correlation length exponent, 
i.e., Tsg/J ~ 0.19 and vsg = 1-2. The best scaling is ob- 
tained with ?7sg = —0.1. 
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FIG. 9: Finite-size-scaling plot of the squared spin-glass order 
parameter with assuming a zero-temperature spin-glass tran- 
sition Tsg = 0. In the main panel, only the data for larger 
sizes and at lower temperatures, i.e., those at T/J < 0.19 and 
with L > 12, are plotted. In the inset, the same scaling plot 
using all the data is shown. 



uses in the scaling plot the data at low temperatures, 
lower than the CG transition temperature T/J < 0.19, 
and the data for larger sizes L > 12, the scaling turns 
out to work well: See the main panel. By contrast, If 
one includes in the scaling plot the data for the small- 
est size L — 8 and at high temperatures T/J > 0.19, 
a significant deviation from the scaling is observed for 
these data: See the inset. In sharp contrast to the scal- 
ing plot of Fig. |S] with Tsg = Tog, we have observed 
here a better scaling for larger sizes, and a systematic 
deviation from the scaling for smaller lattices. In that 
sense, the present finite-size scaling analysis is fully con- 
sistent with the occurrence of a T = SG transition, 



as has long been believed in the community4i£*§*£i2*2*i2i22i 

Furthermore, the exponent associated with the possible 
T = SG transition happens to be rather close to the 
previous estimates based on the numerical domain-wall 
methodic 

Of course, as discussed above, the CG transition occur- 
ring at T/J ~ 0.19 would necessarily affect the nature of 
the SG ordering, even if the Hciscnberg spin itself does 
not order at T = Tcg- Thus, even if the SG transition oc- 
curs only at T = 0, an intrinsic SG critical phenomenon 
associated with this T = transition should set in at 
low temperatures below Tcg, whereas the data at and 
above Tcg would be "contaminated" by the CG transi- 
tion which might well change the ordering behavior of 
the Hciscnberg spin via the associated phase-space nar- 
rowing. 

Hence, although our present data are fully consistent 
with the occurrence of the T = SG transition, in order 
to see such a behavior clearly, one has to choose the scal- 
ing region carefully. Inclusion of the data of smaller sizes 
and at higher temperatures in the analysis would easily 
deteriorate the quality of the scaling plot, leading to the 
opposite conclusion. We believe that this is indeed the 
situation of the recent study of Ref. l39l in which a simul- 
taneous spin and chiral transition at a finite temperature 
Tsg = Tcg was concluded. 



D. Correlation length 

The temperature dependence of the normalized SG and 
CG correlation lengths, £sg/L and £cg/L, for various 
sizes are shown in Fig. ^| In contrast to the Binder 
parameter shown in Fig. [21 the normalized correlation 
lengths of L = 8 and 12 shown in in Fig. II Of a) exhibit 
a clear crossing at a temperature around T/J ~ 0.2 for 
both cases of the SG and the CG. The observed behavior 
is consistent with the behavior recently reported by Lee 
and Young2£ for the 3D Heisenberg SG model with the 
Gaussian coupling for the sizes up to L = 12. 

We now extend the system size up to L = 20, and 
the result is presented in Fig. HOf b) . While the crossing 
temperature for larger sizes L = 16 and 20 shifts to- 
ward lower temperature for both cases of the SG and the 
CG, the CG correlation length still has a clear crossing 
around T/J ~ 0.19, very close to the estimate of Tcg in 
the previous subsections, with a finite crossing angle. On 
the other hand, for the SG correlation length, the cross- 
ing becomes weaker and almost fades away. Namely, the 
curves of L = 16 and 20 merge nearly tangentially with 
a vanishing crossing angle. The L = 16 and 20 curves 
of (,sg/L stay on top of each other in an entire tem- 
perature region studied below T/J ~ 0.19, as if they 
were in the critical KT-like phase. Hence, for the SG, 
with increasing L, not simply the crossing temperature 
shifts toward lower temperature, but the crossing-angle 
becomes smaller and almost vanishes. This is in contrast 
to the behavior of the CG correlation length where the 
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FIG. 10: The correlation length divided by the linear size 
L plotted against the temperature for the chiral-glass £cg / L 
(filled symbols) and for the spin-glass £sg/L (open symbols). 
The data for smaller sizes [L = 8 and 12) are shown in the 
upper figure, and those for larger sizes (L — 16 and 20) are 
shown in the lower figure. 



crossing-angle remains finite with increasing L. It thus 
seems possible that the further increase of L eventually 
leads to the disappearance of the crossing for £sg/L, at 
least in the temperature range studied here. 

This would be consistent with the size-crossover ex- 
pected from the spin-chirality coupling/decoupling pic- 
ture, and with our observation in Sec. IV Bl that the de- 
coupling length scale is about L = 20. Unfortunately, 
at present, we cannot go to lattices larger than L = 20 
due to the limitation of our computation capability. We 
certainly expect, however, that the crossing of £,sg/L 
eventually disappears, or at least shifts to a temperature 
considerably lower than the CG transition temperature 
Tcg/J — 0.19, if we could study lattices considerably 
larger than L = 20. For now, we only mention that, 
although the recent data of the normalized correlation 
length for smaller lattices of L < 12 might look rather 
conclusive at first sightj2£ in view of our present data for 
larger sizes presented in Fig. HOf b). it is still difficult to 
draw a definite conclusion about the ordering nature of 



the model only through the correlation-length measure- 
ments. 



E. Overlap distribution 

In Fig. 1111 we show the chiral-ovcrlap distribution func- 
tion; upper figure (a), and the diagonal-spin-overlap dis- 
tribution function; lower figure (b), at a temperature 
T/J = 0.15 below the CG transition temperature. One 
sees from Fig. Illf a) of the chiral-overlap distribution 
function P(q x ) that, with increasing L, the side peaks 
corresponding to the CG EA order parameter q^Q grow 
and sharpen, which indicates the occurrence of the CG 
long-range order. In addition, a central peak at q x = 
shows up for L > 12, which also grows and sharpens 
with increasing L. The existence of this central peak co- 
existing with the side peaks suggests the occurrence of a 
one-step-like RSB in the CG ordered state. This feature 
is also consistent with the existence of a negative dip in 
the CG Binder parameter gcG an d with the crossing of 
gco occurring on the negative side, as was discussed in 
Sec. IV Al The behavior of P(q x ) observed here is similar 
to the previous reports for the 3D Heisenberg SG with 
the Gaussian coupling^ and the related Heisenberg SG 
modelsi 13 i 22 i 38 By contrast, such a one-step-like feature 
of the overlap distribution has never seen in the Ising SGs 
both with the short-range^i and infinite-range^ interac- 
tions, nor in the Heisenberg SG with the infinite-range 
interaction 2 ^. 

Fig- Hlf b) represents the size dependence of the 
diagonal-spin-overlap distribution function P(gdiag) de- 
fined by Eq. (|18|l . For larger L > 16, the distribution 
function P(qdiag) has only a single peak at <7diag = 0, 
which grows with increasing L, without any other diver- 
gent peak. This is in sharp contrast to the triple-peak 
structure observed in the chiral-overlap distribution func- 
tion P(q x ) of Fig. rnT a). peaked at q x = and ±9cg- ^ 1S 
also in contrast to the double-peak structure observed in 
f (gdiag) of the mean-field Heisenberg SK model, peaked 
at q = ±^qEA~ As discussed in Sec. II VI since the di- 
verging peaks at gdiag = ±j9ea should arise in P((?diag) 
in the possible SG ordered state with a non-zero EA SG 
order parameter^ the absence of any divergent peak at 
non-zero Qdiag for larger L strongly suggests that the sys- 
tem is in the SG disordered state even at this low tem- 
perature T/J = 0.15. Interestingly, a closer inspection 
of Fig. Illf b') reveals that a weak double-peak structure 
can be seen for smaller sizes corresponding to the spin- 
chirality coupling regime, L = 8 and 12. However, such 
a double-peak structure in P(<?diag) tends to disappear 
for larger sizes corresponding to the spin-chirality decou- 
pling regime, L > 16. Again, this could be interpreted as 
the size-crossover from the small-size pseudo SG order to 
the large-size SG disorder, as is naturally expected from 
the spin-chirality coupling/decoupling picture. 



12 



5r- 



0.035 



0.03 - 



0.025 



0.02 



0.015 



0.01 - 



0.005 - 




-0.06 -0.04 -0.02 0.02 0.04 0.06 




FIG. 11: The size dependence of the chiral-overlap distri- 
bution function; upper figure (a), and of the diagonal-spin- 
overlap distribution function; lower figure (b), at the lowest 
temperature of the present simulation, T/J = 0.15. 



F. Equilibrium auto-correlation functions 

Next, we discuss the ordering behavior of the model 
by studying its equilibrium dynamics. In Fig. 1121 we 
show the MC time dependence of the chiral and spin 
autocorrelation functions for our two largest sizes L = 
16 and 20. Here, the time is measured in units of the 
standard heat-bath MC steps without the temperature- 
exchange procedure. 

In the chiral autocorrelation function C x (t) shown in 
Fig. I12f a). no appreciable difference is observed between 
the data of L = 16 and 20 in the time window of t < 
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FIG. 12: The Monte Carlo time dependence of the chiral 
autocorrelation function C x (t); upper figure (a), and of the 
spin autocorrelation function C s {t); lower figure (b), at vari- 
ous temperatures both below and above Tea/ J — 0.19. The 
system size is L = 16 (given by symbols) and L = 20 (given 
by thin lines). The inset is an enlarged view of C 3 (t) in the 
short-time region where the finite-size effect is negligible. 



10 4 , beyond which a weak size effect is appreciable. The 
spin autocorrelation function C s (t), by contrast, is more 
susceptible to the finite-size effect, as can be seen from 
Fig. I12r b). Even in this case, however, the data in the 
time window t < 10 3 shows a negligible size effect as 
shown in the inset. 

As can be seen from Fig. I12f a). C x (t) shows a down- 
ward trend above T/J = 0.19, an upward trend below 
T/J = 0.19, and a near linear behavior at T/J = 0.19. 
In order to quantify this, we fit the data at each tem- 
perature to the form H27|) and plot the x 2 -deviation of 
the fit in Fig. ^| as a function of the temperature around 
the expected CG transition temperature. The plot has a 
minimum around T/J = 0.19(1), at which the data are 
optimally fitted to a power law. This estimate of Tqg 
based on the chiral autocorrelation function agrees well 
with those obtained from the CG Binder parameter and 
the CG order parameter. 

We also test a dynamical scaling analysis of the chi- 
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FIG. 13: Residuals per degrees of freedom associated with 
the x 2 -fitting of the chiral autocorrelation function (marked 
by filled circle), and an estimated effective exponent P/zv 
(marked by open circle) plotted against the temperature. 



FIG. 14: Dynamical scaling plot of the chiral autocorrela- 
tion function. The best scaling is obtained with choosing 
Tcg = 0.19, Pcg = 0.8 and zcg^cg = 5.4. The upper branch 
represents the scaling of the data below Tcg, while the lower 
branch represents that above Tcg- 



ral autocorrelation function C x (t). As can be seen from 
Fig. El the dynamical scaling works well both above and 
below Tcg, with the scaling parameters Tcg/ J = 0.195, 
Pcg = 0.8 and zcg^cg = 5-0- The present estimate of 
@cg is slightly smaller than, but is not far from the pre- 
vious estimate of Ref. f° r the 3D Heiscnberg SG with 
the Gaussian coupling (3cg — 1-1. 

As can be seen from Fig. I12f b). by contrast, the 
spin autocorrelation O s (t) shows a downward trend at 
longer times at any temperature studied, suggesting an 
exponential-like decay characteristic of the disordered 
phase. A closer inspection of the data of C s (t) reveals 
that the data below Tcg exhibit a weak hump-like struc- 
ture at short times t ~ 10 2 , though this hump eventually 
gives way to the down-bending trend characteristic of the 
disordered phase at longer times t > 10 3 : See the inset of 
Fig. I12f b). This hump- like structure observed in C s (t) 
at short times might be a manifestation of the trivial 
spin-chirality coupling expected at short time scales, and 
is not likely to be an indication of the SG long-range or- 
der, since the downward trend is recovered at longer time 
scales: See Sec. I VII for further details. Hence, from our 
dynamical data, we conclude again that the CG transi- 
tion occurs at Tcg/ J = 0.19(1), without accompanying 
the simultaneous SG order. 



G. Equilibrium correlation time 

Generally, the temporal decay of the autocorrelation is 
characterized by the temperature-dependent characteris- 
tic time scale, the correlation time t(T), which represents 
a dynamical crossover from the short-time critical behav- 
ior to the long-time relaxation. The correlation time t(T) 
diverges as r(T) ~ |T — T c \~ zv when the temperature 
approaches T c from above. One promising method of es- 
timating t from the autocorrelation function has been 
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FIG. 15: Scaling plot of the chiral ratio function; upper figure 
(a), and of the spin ratio function; lower figure (b), at various 
temperatures for the sizes L — 16 and 20. 
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proposed by Bhatt and Young^, who employed a scaling 
analysis of the dynamical ratio function. For the spin 
autocorrelation, this reads as 



R a (t) 



C s (t) 



(29) 



The corresponding chiral ratio function is defined by 

o x (t) 



R x {t) 



(30) 



Because the ratio function is dimensionless, the prefactor 
t-P/zv m Eq, ||28|l is canceled out. The dynamical scaling 
form of R(t) is then given as a single- variable function of 
t/r, 



R(t) = R(t/r), 



(31) 



where R is a scaling function. If one appropriately 
chooses the scaling parameter r which depends on the 
temperature and the system size, the ratio functions 
should be scaled on to a single curve. Using this method, 
Bhatt and Young^ successfully estimated the correlation 
time of the short-range Ising EA model and the mean- 
field Ising SK model. Subsequently, this method has been 
extended to non-equilibrium relaxation, where the ratio 
function depends not only on the measurement time t 
but also on the waiting time t w . The off-equilibrium 
method was applied recently by Matsumoto, Hukushima 
and Takayama to the 3D ± J Heisenberg SG4^ 

Here we use this method to estimate the correlation 
times both for the spin and for the chirality. In compari- 
son with the previous off-equilibrium study^ the present 
equilibrium study has an advantage that one needs not 
extrapolate to an equilibrium limit, i.e., needs not take 
the t w — > oo limit. In Fig. ED we show the scaling plot of 
the chiral and spin ratio functions. We note that both the 
spin and chiral scaling functions are described roughly by 
an exponential form. 

In order to compare the spin and chiral correlation 
times, denoted by tsg and tcgj respectively, we plot 
them in Fig. 1161 as a function of the temperature. In 
the figure, we combine the data with those obtained 
in a wider temperature range by off-equilibrium simu- 
lation of Ref. |4J. As can be seen from Fig. the 
chiral correlation time is shorter than the spin correla- 
tion time at higher temperatures, similarly to the be- 
havior of the correlation length discussed in Sec. IV Dl 
As the temperature is decreased, the chiral correlation 
time tcg grows faster than the spin correlation time tsg, 
and eventually exceeds t§q at a certain characteristic 
temperature T = T x (L) . The size dependence of this 
crossover temperature T = T x (L) is apparently weak: 
We get T = T x ~ 0.24 both for L = 16 and 20. It 
strongly suggests that, even in the thermodynamic limit, 
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FIG. 16: The temperature dependence of the the chiral and 
spin correlation times for the sizes L — 16 and 20. The cor- 
responding data obtained from the off-equilibrium simulation 
of Ref. 113 are also included. 



the chiral correlation time exceeds the spin correlation 
time at a crossover temperature T x (oo), which is lo- 
cated somewhat above the CG transition temperature 
Tcg/ J — 0.19. It means that, with decreasing the tem- 
perature, the relevant degree of freedom dominating the 
long-time ordering behavior changes from the spin to the 
chirality at T = T x . In order to further illustrate this 
changeover, we show in Fig. 1171 the time-dependence of 
the spin and chiral ratio functions at two representative 
temperatures T/J = 0.25 and 0.20, each above and below 
T x . As can be clearly seen from the figure, with decreas- 
ing the temperature across T x , the temporal decay of the 
chiral ratio function becomes much slower than that of 
the spin ratio function: Compare the two arrows in the 
figure. 

The time scale associated with such a crossover, £ x , is 
roughly estimated to be 10 5 ~ 10 6 MCS. For more precise 
estimate of i x , more quantitative analysis of the size de- 
pendence of the crossover time scale would be necessary. 
Naturally, this crossover time t x gives a measure of the 
time scale above which the spin-chirality decoupling can 
be observed in dynamics. Thus, the spin-chirality decou- 
pling in the dynamics would be eminent only at temper- 
atures lower than T x /J ~ 0.24 and at times longer than 
t x ~ 10 5 — 10 6 MCS. This crossover time scale is rather 
long, yet, is finite. It is important to realize that, in order 
to resolve the controversy concerning the presence or the 
absence of the spin-chirality decoupling in the Heisenberg 
SG, one has to probe the equilibrium dynamics beyond 
this crossover time scale 10 5 — 10 6 MCS at temperatures 
lower than T x /J ~ 0.24, about some 20 ~ 30 % above 

rp 45 

As argued in Sec. ^ from the spin-chirality cou- 
pling/decoupling picture, a similar phenomenon is ex- 
pected also in the spatial correlation of the model in 
terms of the length scale. Namely, one expects that 
at a certain crossover temperature T' x , which is prob- 
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FIG. 17: Temporal decay of the chiral and spin ratio func- 
tions at temperatures T/J — 0.25 and 0.20. The system size 
is L = 16. With decreasing the temperature from T/J — 0.25 
to 0.20, the chiral relaxation slows down much more slowly 
than the spin relaxation, as illustrated by the arrows. 



FIG. 18: The temperature dependence of the chiral and spin 
correlation lengths for finite systems. The system size is L = 
16 and 20. The data are the same as those shown in Fig. 1101 
but not divided by L here. The curves are polynomial fits 
of the data which are extrapolated to lower temperatures to 
deduce the crossing temperature given in the text. 



ably close to the dynamical crossover temperature T x 
discussed above, the CG correlation length £cg exceeds 
the SG correlation length £sg- This changeover of the 
dominant length scale gives a crossover length scale L x 
above which the spin-chirality decoupling is eminent in 
spatial correlations. Unfortunately, unlike the case of 
the correlation time, the limitation of the available sys- 
tem size prevents us from directly estimating L x . In 
Fig. El we plot the temperature dependence of the CG 
and SG correlation lengths for the sizes L = 16 and 20. 
For these sizes, the crossing of £sg and £cg occurs at a 
temperature lower than the CG transition temperature, 
in contrast to the case of the correlation time. Neverthe- 
less, the crossover temperature at which £sg and £cg of 
finite L cross, tends to increase with increasing L. If we 
roughly estimate the crossover length scale of finite sys- 
tems by extrapolating the data of Fig. El we tentatively 
get L x ~ 11 (L = 16) and L x ~ 14 (L = 20). These 
results are certainly not inconsistent with our estimate of 
L x ~ 20 based on the behaviors of the SG order param- 
eter, the dimcnsionless correlation length £/L and other 
quantities. 

H. The A and G parameters 

We have also calculated the A and G parameters de- 
fined in Sec. IIVI both for the CG and SG orders. In 
Fig. El the temperature and size dependence of the A 
and G parameters for th CG order, Aqg and Gqg, is 
shown. Although the data are rather noisy due to the 
large sample-to-sample fluctuations, the A parameter of 
different L show a crossing and a peak around the ex- 
pected CG transition point T/J ~ 0.19, as can be seen 
from Fig. I19f a). In particular, with increasing L, Aqg 
stays non-zero below Tqg, indicating that the CG or- 
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FIG. 19: The temperature and size dependence of the A 
parameter of the chirality; upper figure (a), and of the G 
parameter of the chirality; lower figure (b). 
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FIG. 20: The temperature and size dependence of the A pa- 
rameter of the spin; upper figure (a), and of the G parameter 
of the spin; lower figure (b). 



dered state is non-self-averaging. These findings com- 
bined with the peculiar shape of P(q x ) shown in Sec. IV El 
suggest that the CG ordered phase accompanies an RSB 
with a non-self-averagencss. For the corresponding G 
parameter, the crossing is not so clear, as is shown in 

Fig. cab). 

The temperature and size dependence of the A and G 
parameters for th SG order, A^g and Gsg ; is shown in 
Fig. [20| As shown in Fig. ^Of a - ). the A parameter ex- 
hibits a crossing, although the crossing temperature is 
located considerably above T c for this range of L. One 
might be tempted to interpret such a crossing of A$q as 
an unambiguous evidence of the occurrence of the stan- 
dard SG transition. However, one has to be careful here: 
Although the crossing of A$g is certainly a signature of 
some sort of phase transition occurring there, it does not 
necessarily mean the occurrence of the standard SG tran- 
sition characterized by a non-zero SG order parameter. 
A non-zero Aqq persisting in the L — > oo limit simply 
means that the SG order parameter q$Q, or the SG sus- 

(2) 

ccptibility xsg = Nq^Q, is non-self-averaging. Below 
the CG transition temperature, one expects that the SG 



order parameter is still Gaussian-distributed around zero 
with a width corresponding to the finite SG susceptibility 
Xsg j while the width exhibits sample-to-sample fluctua- 
tions leading to the non-self-averaging xsg- The latter 
is a natural consequence of the phase-space narrowing 
which should inevitably accompany the CG transition 
with the one-step-likc RSB. Hence, the crossing of Age 
and a finite ^4sg remaining in the L — > oo limit below 
Tec are compatible with the absence of the standard 
SG long-range order, and is entirely consistent with the 
the CG transition not accompanying the standard SG 
long-range order. 

As shown in Fig. I20f b~) , the G parameter of the spin 
exhibits a crossing around 7c g- The relation Eq. J23J, 
combined with our observation in Fig. |3fb), indicates 
that Gsg also takes a non-zero value below Tcg- Thus, 
the observed crossing of Gsg is just as one expects for 
the CG transition. In other words, one cannot interpret 
the crossing of Gsg a s an indicator of the onset of the 
standard SG long-range order. 



VI. DISCUSSIONS 

In this section, in view of our MC results presented in 
the previous section, we wish to examine and discuss the 
recent numerical studies on the 3D Heisenberg SG. Many 
of these studies suggested, contrary to our present study, 
that the spin and the chirality ordered simultaneously at 
a finite temperature with a common correlation length 
exponent z^sg = ^cg> *- e -, no spin-chirality decoupling 
in the 3D Heisenberg SGi 16 i 18 i 19 i 20 i 39 i 46 Below, we wish 
to make some comments on these numerical works from 
the standpoint of the spin-chirality coupling/decoupling 
picture. 



A. Stiffness method 

First, we wish to discuss the analyses based on the 
stiffness methodiiSiiLil In this method, one computes 
by some numerical means the change of the ground-state 
energy of finite systems of size L under the appropriate 
change of boundary conditions imposed on the system. 
This energy is called a stiffness energy (or a domain-wall 
energy), AEl, which gives a measure of an energy scale 
of low-energy excitations of size L. For large L, AEl is 
expected to behave as a power- law, AEl ~ L e , 9 being 
a stiffness exponent. If 9 < 0, the system remains in the 
disordered state at any nonzero temperature, whereas if 
9 > the system possesses a finite long-range order at 
low enough temperatures with T c > 0. Here, we discuss 
this method first in conjunction with the detection of the 
standard SG order, leaving the detection of the CG order 
later. 

The nontrivial part of this stiffness method concerns 
with the choice of the boundary conditions employed in 
computing the stiffness energy. There could be various 
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choices, and the behavior of AEl might in principle de- 
pend on these choices particularly for small sizes accessi- 
ble in numerical simulations. The most standard choice 
is the combination of the periodic and the antipcriodic 
boundary conditions (P/AP). In the case of the Heisen- 
berg SG, the P/AP combination necessarily accompanies 
a flipping of the chirality (remember that the chirality 
of the Heisenberg spin is odd under the spin inversion 
S — > —S), so that the P/AP combination should detect 
the chiral order for large enough In order to detect 
the standard SG order independently of the CG order 
by this stiffness method, Ref. introduced the "rota- 
tion" boundary conditions (ROT), which imposed a 7r 
rotation on the boundary spins without flipping the chi- 
rality, which was combined with the standard P bound- 
ary conditions in calculating the stiffness energy^ Such 
a P/ROT combination applied to the 3D Heisenberg SG 
yielded a negative 8, i.e., 8 ~ —0.51 for the Gaussian 
coupling, and 8 ~ —0.49 for the ± J coupling, which im- 
plied the absence of the standard SG order at nonzero 
temperature^ 

By contrast, Matsubara, Endoh and Shirakura pro- 
posed a different choice of boundary conditions in com- 
puting AEl, i.e., to use the free (open) boundary con- 
ditions as a reference and impose the rotational-twist 
to such "optimized" spin configurations obtained under 
the free boundary conditions in which the stress at the 
boundary is released i&lL These authors observed that 
the stiffness exponent evaluated in this way was largely 
positive, close to the spin- wave exponent 8 = 1, and 
argued that the 3D Heisenberg SG exhibited a finite- 
temperature SG transition. The method similar in spirit 
to the one used in Refs. IT^ andll7l was also applied to the 
XY SG by Kosterlitz and Akino^i leading to the similar 
conclusion. Thus, the result obtained by applying the 
free/twisted- free (F/TF) boundary conditions, 8 > im- 
plying Tsg > 0, is in sharp contrast to the result obtained 
by applying the P/ROT boundary conditions, 8 < im- 
plying Tsg = 0. Discrepancy between the stiffness ex- 
ponents evaluated by the different choices of boundary 
conditions was also observed in other SG models, e.g., in 
the 2D Ising SGi£ 

Although the authors of Rcfs. [Hand [13 argued that 
their "optimized" boundary conditions were superior to 
the other choices, a theoretical basis of such a claim seems 
not so obvious. For example, one might make the follow- 
ing counter-argument that one should not optimize the 
boundary conditions in calculating AEl'- In the spirit of 
the domain- wall renormalization-group idea by Bray and 
Moored, AE(L) represents an energy scale associated 
with the interaction between the coarse-grained blocks 
of size L in an infinite SG sample. Since these blocks are 
necessarily subject to the strong frustration effect caused 
by the interaction with the neighboring blocks surround- 
ing them, an optimization of their energy, independently 
at each block ignoring the inevitable frustration effect due 
to the interaction with the neighboring blocks, is hardly 
compatible with the original RG idea. One may thus 



argue that, in calculating AEl in SGs, one should not 
make the optimization of boundary conditions referring 
to the particular bond realization of each sample. 

Concerning the apparent difference of the stiffness ex- 
ponents arising from the different choices of boundary 
conditions, there generally exist two possibilities: Either, 
(i) the observed difference is a finite-size effect where 
there is a single stiffness exponent for large enough lat- 
tices independent of applied boundary conditions, or, (ii) 
the observed difference is a bulk effect which persists even 
in large enough lattices. In the case of the 2D Ising SG, 
Carter, Bray and Moore numerically observed that, al- 
though both the P/AP and F/TF boundary conditions 
yielded the same stiffness exponent asymptotically for 
large L, i.e., the possibility (i) above, the finite-size ef- 
fect was much reduced in the P/AP than in the F/TF^S 
For vector SG, there so far exists no convincing evidence 
which of the above (i) and (ii) is really the case. In 
any case, a practical question we face with is which set 
of boundary conditions gives a true asymptotic answer 
from smaller sizes accessible in simulations. 

One plausible criterion might be that, among all pos- 
sible excitations in the system, the one giving the low- 
est excitation energy AEl, or equivalcntly, the one giv- 
ing the smallest stiffness exponent 8, should be chosen. 
The reason is simply because among all possible excita- 
tions the one with the lowest excitation energy should 
be most efficient in destroying the order so long as it 
has non-negligible weight in the thermodynamics, and 
would dominate the low-energy dynamics of the model. 
Under this criterion, when the different choices of bound- 
ary conditions yield different 8 values, the one giving the 
smallest 8, or the most negative 8, should be chosen. In 
particular, when one set of boundary conditions yields 
a positive 8 while the other yields a negative 8, the one 
giving a negative 8 should be chosen. If so, in the case of 
the 3D Heisenberg SG of our interest, the P/ROT com- 
bination without any optimization procedure should be 
chosen since it gives the lowest 8 (negative 8) reported 
so fariS This suggests that the standard SG order in the 
3D Heisenberg SG occurs only at T = 0. Of course, it 
is still possible that some other type of boundary condi- 
tions might yield a still smaller 8, but it does not change 
the conclusion that the SG order occurs only at T = 0. 

It should also be remembered that the types of low- 
energy excitations generated via a particular choice of 
boundary conditions is only a subset of all possible exci- 
tations in the system: They are basically wall-like exci- 
tations, not including more complex excitations like, say, 
a "vortex" excitation which is possible in the Heisenberg 
SG reflecting the SO(3) nature of its order parameter 
spaced or a "sponge" excitation which is closely related 
to the RSB structure of the ordered state^S Unfortu- 
nately, we have little knowledge concerning what is the 
most relevant low-energy excitation governing the order- 
ing of the system, and hence, have no well-based criterion 
to choose one set of boundary conditions from the others 
as superior. Although we feel that our argument above 
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speaks for a zcro-tcmperature SG transition in the 3D 
Heisenberg SG, it would be fair to say at present that no 
definitive conclusion can be drawn solely based on this 
stiffness method. 

We finally wish to refer to the stiffness method in de- 
tecting the chiral order. As mentioned, since the sign of 
the chirality is flipped by changing the boundary condi- 
tions from P to AP, the most standard P/AP combina- 
tion could be used in detecting the chiral order, at least 
for large enough L. In practice, however, the applica- 
tion of the reflecting (R) boundary is more efficient in 
detecting the chiral order, as shown in Ref. 0. The chiral 
stiffness exponent of the 3D Heisenberg SG determined in 
this way turned out to be positive implying a CG transi- 
tion occurring at a nonzero temperature £ Other aut hors 
also reported a positive value for the chiral stiffness ex- 
ponent both for the 3D XY SG~ and the 3D Heisenberg 
SGil. 



such a hump of C s {t) appears only at times shorter 
than the crossover time scale t x (t x was estimated to 
be t x « 10 5 - 10 6 above T CG ), while C s (t) eventually de- 
cays at long enough times t >> t x Indeed, as was shown 
in the inset of Fig. 1121 such a hump-like weak structure 
of the spin autocorrelation was discernible in our present 
data of C s (t) at short times t ~ 10 2 , which, however, 
eventually decayed toward zero at longer times. 

Bcrthier and Young also observed in their recent off- 
equilibrium simulation of the 3D Heisenberg SG a weak 
hump-like structure in the spin autocorrelation in the 
time range t < 10 4 , which corresponded to the quasi- 
cquilibrium regime 4 ^. These authors interpreted the ob- 
served hump as an evidence of a nonzero SG long-range 
order at that temperature. As noted above, however, 
such a hump is also consistent with the the spin-chirality 
coupling/decoupling picture as long as the hump is ob- 
served only at shorter times t < t x . 



B. Equilibrium dynamics 

Matsubara, Shirakura and Endoh reported a further 
evidence of the simultaneous spin and chiral transition in 
the 3D Heisenberg SG by investigating the equilibrium 
spin dynamics^. In order to eliminate the effect of global 
spin rotations inherent to finite systems, Matsubara et al 
introduced an artificial global-rotation correction in the 
spin dynamics of the model. They observed that the 
modified spin autocorrelations adjusted by the global- 
rotation correction exhibited at lower temperatures a ten- 
dency to approach a nonzero value at longer times, which 
was interpreted as an evidence of a finite SG long-range 
order. 

It should be noticed here that, when one looks at a 
quantity which is even under the symmetry transforma- 
tion of the Hamiltonian like the modified spin autocorre- 
lation function of Ref . ITsl one needs to examine its size 
dependence carefully. As is well-known, an even quan- 
tity in finite systems always takes a nonzero value even 
above T c due to the finite-size effect, where this nonzero 
value decreases with the size L, eventually vanishing as 
L — > oo above T c . (Indeed, in an extreme occasion of 
a single spin, the modified spin autocorrelation function 
as computed by Matsubara et al does not decay at all 
even at an infinite temperature !) The ordering behavior 
of the modified spin autocorrelation as observed by Mat- 
subara et al might possibly be caused by the finite-size 
effect. In order to refute such suspicion, one needs to 
study its size dependence carefully, whereas the analysis 
of Ref. ITsl was limited to a fixed size L — 16. 

It should be stressed that, even within the spin- 
chirality coupling/decoupling scheme, it is still possible 
that the spin autocorrelation function C s (t) of an infinite 
system exhibits below Tqg & hump-like weak structure 
at short times as illustrated in Fig[3 which is an echo 
of the plateau-like structure of the chiral autocorrelation 
function. In the temperature range Tsg < T < Tec 



C. Nonequilibrium dynamics 

Nakamura and Endoh applied a non-equilibrium 
method to study the SG and the CG orderings of the 3D 
± J Heisenberg SG>i^ Analyzing the time dependence of 
the initial growth of the SG and the CG susceptibilities 
with use of a dynamical scaling, these authors concluded 
that the spin and the chirality ordered simultaneously 
at a finite temperature T/J = 0.21 ~ 0.22. While the 
lattice size studied L < 59 was rather large, the crucial 
question to be addressed is whether the long-time limit 
t — » oo was safely taken justifying the use of the dynami- 
cal scaling. In other words, although the nominal lattice 
size studied was large, the equilibrated length scale ac- 
tually probed in these off-equilibrium simulations might 
be rather short. In fact, their non-equilibrium method 
is uncontrolled with respect to the time scale toward 
equilibrium. Since the equilibration time could easily 
become a huge number in SG, care has to be taken as 
regards the equilibrated length scale actually probed by 
this type of non-equilibrium simulation. As one judges 
from the maximum values of the SG and the CG sus- 
ceptibilities reached by the off-equilibrium simulation of 
Ref.ll9l the "dynamical correlation length" still remained 
rather short: Namely, even around the transition tem- 
perature Tqg/J — 0.2, it reached around 10 lattice spac- 
ings for the spin, and only 1 or 2 lattice spacings for 
the chirality. The dynamical chiral correlation length 
stayed particularly short. This is consistent with a re- 
cent off-equilibrium simulation by Berthier and Young 
in which the dynamical chiral correlation length stayed 
much shorter than the dynamical spin correlation length 
in the investigated time ranged. Here note that, irre- 
spective of the question of whether the CG transition 
accompanies the simultaneous SG transition or not, the 
chiral correlation length in equilibrium should diverge at 
and below the CG transition temperature Tcg/J — 0.2. 
Hence, the observation above simply tells that, even at 
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the maximum simulation time of Refs^] andl46l the sys- 
tem still stayed in an extreme initial time regime. In or- 
der to deduce the equilibrium ordering properties from 
these off-equilibrium data, one is forced to extrapolate 
the behavior around £cg ~ 1 to £cg = 00 7 which could 
be dangerous in the present model since the model might 
possess the characteristic crossover length scale at around 
20. 

One may feel that the dynamical spin correlation 
length reached in the off-equilibrium simulation of 
Ref. E £ » 10, might be reasonably large for deduc- 
ing the ordering properties of the spin. However, we feel 
it is not enough. This length scale of 10 is still not large 
enough compared with the crossover length scale esti- 
mated in the present work L x « 20. Remember that the 
spin-chirality decoupling, if any, is a long scale phenom- 
ena observable at length scale longer than L x . Second, 
in the off-equilibrium simulations of Refs.ll9land l46l even 
when the dynamical SG correlation length grows around 
10 lattice spacings, the Zi chiral degree of freedom was 
not equilibrated at all at this length scale, in sharp con- 
trast to the fully equilibrated simulation as was done in 
the present paper. In other words, at the length scale of 
10 lattice spacings, the chiralities arc little thermalizcd 
and are virtually frozen in the non-equilibrium pattern, 
while only the SG correlation grows modestly in such a 
noncquilibrium chiral environment. After all, however, 
we have to understand the spin dynamics at long enough 
length scales at which the chiral degree of freedom is 
also fully thermalizcd. Thus, the spin dynamics as ob- 
served in the off-equilibrium situations of Refs. and 
l46l may not faithfully represent the close-to-equilibrium 
critical dynamics of the original model. 

A similar dynamical simulation on the 3D ± J Heisen- 
berg SG was performed by Matsumoto, Hukushima and 
Takayama44 They also made a dynamical scaling analy- 
sis taking the effect of global spin rotations into account. 
In this study, the time scale toward equilibrium was con- 
trolled via the analysis of the wai ting -time dependence of 
the results. In contrast to Refs. Il9l and liti Matsumoto 
et al suggested that their data were consistent with the 
separate spin and chiral transition, i.e., Tsg < ?cg- 

Berthier and Young also argued that their observa- 
tion that the dynamical CG correlation length stayed 
shorter than the dynamical SG correlation length pre- 
sented an evidence that the spin, rather than the chiral- 
ity, was the order parameter of the transition^ Some 
caution is required in drawing the final conclusion from 
this observation, though. Both the spin and the chi- 
rality length scales pro bed by the off-equilibrium simula- 
tions of Refs^^ andl46l are still shorter than the crossover 
length estimated in the present paper £ x — 20. Hence, 
within the spin-chirality coupling/decoupling scenario, 
there still exists a good possibility that the dynamical 
CG correlation length eventually outgrows the dynami- 
cal SG correlation length at longer times. 

In Ref. l46l the aging phenomena were persistently ob- 
served at lower temperatures, not only for the chirality 



but also for the spin, which was interpreted as an evi- 
dence of the occurrence of simultaneous spin and chiral 
transition^ Again, this cannot be taken as an unam- 
biguous indicator of a finite-temperature SG transition, 
since the aging phenomena could arise simply when the 
time scale of measurements becomes comparable to the 
longest relaxation time in the system which could be ex- 
tremely long in SG even in the paramagnetic phase. For 
example, in the 2D Ising SG which is known to exhibit no 
finite-temperature SG transition, clear aging phenomena 
have been observed both in numerical simulations^^ 
and in experiments^^. 



D. Correlation length 

In Ref. [53, Lee and Young calculated by means of equi- 
librium simulations both the SG and the CG correlation 
lengths of the 3D Hciscnbcrg SG with the Gaussian cou- 
pling in the range of sizes 4 < L < 12. Lee and Young ob- 
served a crossing of the dimensionless correlation lengths 
t;/L for different L for both cases of the spin and the chi- 
rality, and concluded that the spin and the chirality or- 
dered simultaneously at a finite temperature T / J ~ 0.15. 
The behavior of £/L observed in Ref. [2(] turned out to 
be quite different from that of some other dimension- 
less quantities, e.g., the Binder ratio, whereas Lee and 
Young argued that the correlation length might be the 
most trustable quantity to look at. Generally speaking, 
however, £/L is also subject to significant finite-size ef- 
fects, sometimes no better than other quantities^ 

We note that the numerical data of Ref. |2(] are basi- 
cally consistent with our present data for smaller sizes 
L < 12 : See Fig. HUa). As emphasized in Sec. IVDl 
of our present paper, however, the crossing behavior of 
the dimensionless SG correlation length £sg/£ tends to 
change for larger lattices L > 12: The crossing becomes 
weaker and weaker, and £,sg/L of L = 16 and that of 
L — 20 do not quite cross with a finite crossing angle as 
occurs for smaller lattices L < 12, but instead, merge al- 
most tangcntially and stay on top of each other at lower 
temperatures: See Fig. HUf b). In contrast, the dimen- 
sionless CG correlation length £cg/L of L = 16 and that 
of L = 20 persistently exhibits a clear crossing. If such 
a tendency continues for larger lattices, the crossing of 
£,sg/L might no longer occur for large enough lattices, 
at least at the crossing temperature of £cg/L. We thus 
suspect that the crossing behavior of £sg /L as reported 
in Ref. [2(] might be a transient behavior due to the small 
sizes, which reflected the "coupling" behavior expected 
at L < L x w 20. Namely, within the spin-chirality cou- 
pling/decoupling scheme, the SG correlation exhibits the 
ordering behavior similar to the CG correlation at shorter 
length scales at which the spin is coupled to the chiral- 
ity, while at longer length scales at which the spin is 
decoupled from the chirality, the SG correlation eventu- 
ally exhibits the non-ordering behavior different from the 
CG correlation. Unfortunately, the largest lattice size ac- 
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ccssible by the present computational capability L ~ 20, 
being only comparable to the crossover length for the 
spin-chirality coupling/decoupling to occur, is still not 
large enough to clear see this behavior. We do expect, 
however, that the correlation lengths for larger lattices 
L > 20 would eventually exhibit a clear spin-chirality 
decoupling behavior. 

We note that such a coupling/decoupling behavior of 
the SG correlation length in smaller/larger lattices was 
indeed observed recently in the 2D Heisenberg SGi^ For 
the 2D Heisenberg SG with the Gaussian coupling, Kawa- 
mura and Yonchara calculated the dimensionless SG cor- 
relation length £,sg/L up to the size £=40, and found 
that £,sg/L for the smaller sizes L = 10,16,20 crossed 
almost at a common temperature T/J ~ 0.022, dis- 
guising the occurrence of a finitc-tcmpcraturc SG transi- 
tion: See the inset of Fig. 6(a) of Ref. I24I while the data 
for the larger sizes L=30 and 40 data eventually came 
down, no longer making a crossing at T/J ~ 0.022. The 
asymptotic non-ordering behavior observed for L > 20 is 
consistent with a zero-temperature SG transition, which 
has been well established in 2Di2iS4 Meanwhile, the CG 
correlation length exceeds the SG correlation length at 
around T/J ~ 0.022, which might naturally explain the 
reason why £,sg/L for smaller sizes L < 20 exhibited a 
crossing behavior. Anyway, this observation in 2D gives 
us a warning that one should be careful in interpreting 
the crossing-like behavior of £/L observed for smaller 
sizes as an unambiguous evidence of a true SG phase 
transition. 



E. Finite-size scaling of the order parameter 

Matsubara, Shirakura, Endoh and Takahashi made a 
finite-size scaling analysis of the SG order parameter q^Q 
for the 3D ± J Heisenberg SG, and claimed that the qual- 
ity of the scaling was much better when ones assumed a 
nonzero SG transition temperature T$q/J = 0.18 than 
a zero SG transition temperature Tsg/J = Their 
conclusion is in apparent contrast to that of our present 
work based on a similar scaling analysis in Sec. IV Gl We 
note that the quality of the finite-size scaling is some- 
times sensitive to the range of lattice sizes and the range 
of temperatures used in the fit. 

As already noticed in Sec. IV Gl this point could be 
particularly serious in the present model. In the spin- 
chirality coupling/decoupling scheme, the SG correlation 
and the CG correlation are trivially coupled at shorter 
length scale L < L x w 20, so that even the SG order 

(2) 

parameter qg,A would be scaled for smaller sizes with as- 
suming a simultaneous SG and CG transition, with ap- 
parent (not true) SG pseudo-exponents u^q ~ ^cg an d 
1 + ^SG — (1 + ?7cg)/3- Indeed, as was shown in Fig. [HI 
our present data, particularly those of L < 16, turned out 
to be scaled reasonably well by assuming a simultaneous 
SG and CG transition at T/J = 0.19, which we inter- 
preted as a pre-asymptotic pseudo-critical behavior real- 



ized in the short-scale coupling regime. Furthermore, the 
relation between t/cg and ?7gQ mentioned above roughly 
holds at short length scale; (1 + ?7cg)/3 ~ 0.60 versus 
1 + ?7§| ~ 0.88 at T/J = 0.19. At longer length scales, 
however, the spin is eventually decoupled from the chiral- 
ity. Then, if ones continues to put T$g = Tqg m the fit 
of <7gQ, the good data collapse obtained for smaller sizes 
would eventually deteriorate for larger sizes. Indeed, as 
shown in Fig. [8] our L = 20 data of q^l showed such a 
deviation expected for larger sizes. 

Within the spin-chirality coupling/decoupling scheme, 
in order to see the true asymptotic critical behavior of the 
SG transition occurring at Tsg(< Tqq), one has to en- 
ter into the long-scale decoupling regime and well below 
the CG transition temperature, i.e., L > L x m 20 and 
T < Tqg ~ 0.2 J. It should be noticed that, in their scal- 

(2) 

ing fit of <7gQ, Matsubara et al included the data points 
for smaller sizes L = 5, 7, 9 etc, which are expected to 
lie in the short-scale coupling regime, as well as the data 
points at temperatures above Tcg which might lie out- 
side the asymptotic critical regime of the SG transition. 
Hence, the poor scaling reported by Matsubara et al with 
assuming Tqg = might simply be due to the fact that 
the data points used in the fit are not in the correct 
asymptotic regime. 

By contrast, we have observed that, if we use the data 
points only of larger lattices L > 16 and only at tem- 
peratures below T"cg j the data were scaled reasonably 
well even with assuming Tsg = 0: See Fig.|5] Therefore, 
we believe that there still exists a good possibility that 
the SG order occurs only at T = as has widely been 
believed in the community, although it is also quite pos- 
sible that it occurs at a low but nonzero temperature, 
< T SG < T CG - 

As discussed in some detail above, any of the re- 
cent works claiming the simultaneous spin and chiral 
transition in the 3D Heisenberg SG appears not con- 
clusive. As far as the authors are aware, all of these 
observations are consistent with the spin-chirality cou- 
pling/decoupling scheme with the crossover length scale 
of 20 lattice spacings and the crossover time scale of 
10 5 - 10 6 MCS. Rather, we believe that some of the ob- 
servations reported in the present paper give a strong 
numerical support that the SG transition indeed occurs 
at a temperature below the CG transition temperature, 
i.e., T S G < T C g- 

VII. SUMMARY 

In summary, we studied the equilibrium properties 
of the three-dimensional isotropic Heisenberg spin glass 
by means of extensive MC simulations. We presented 
evidence of a finite-temperature CG transition with- 
out accompanying the conventional SG order through 
the observation of various physical quantities including 
the order parameters, equilibrium static and dynamic 
correlation functions, Binder parameters and overlap- 



21 



distribution functions, etc. Our conclusion is in con- 
trast to some of the recent numerical studies on the same 
model, which claimed the simultaneous SG and CG tran- 
sition. We have pointed out that the crossover length 
scale and the crossover time scale associated with the 
spin-chirality coupling/decoupling are crucially impor- 
tant in properly interpreting the numerical data. Around 
the CG transition temperature, these length and time 
scales are roughly estimated to be 20 lattice spacings 
and 10 5 ~ 10 6 MCS, respectively. Below these length 
and time scales, the spin is trivially coupled to the chi- 
rality so that the spin-chirality decoupling, i.e., the SG 
disorder, is difficult to observe. This might give a natu- 
ral interpretation of the discrepancy between our present 
result and the observation of the simultaneous SG and 
CG transition by some other authors. 

Rather, it appears to the authors that our present 
data for larger L are hard to understand based on the 
standpoint of the simultaneous spin and chiral transi- 
tion. Hence, while simulations on still larger lattices with 
L > 20 are required to settle the issue, our present data 



give some support to the spin-chirality decoupling sce- 
nario for the 3D isotropic Heisenberg spin glass. 
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